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Charge induced energy fluctuations in thin organic films: eff^ect on charge transport 

Sergey V. NovikovQ 
^.A'". Frumkin Institute of Electrochemistry, Moscow, Russia 

Effect of static diarges on charge carrier transport in disordered organic materials is considered. 
Long range nature of Coulomb interaction requires to take into consideration a finite thickness of 
the transport layer. Presence of conducting electrodes significantly modifies properties of organic 
medium, removes a long range Coulomb divergence, and makes it possible to calculate in finite form 
statistical properties of organic medium (with nonzero total charge density), relevant for transport 
characteristics. A special attention is paid to the particular case of charge induced disorder - a 
disorder originated from the surface charge located at the rough surface of electrode. We present 
also a generalization of ID model of charge carrier transport to the case of inhomogeneous energetic 
disorder that realizes in for charge induced disorder. 
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I. INTRODUCTION 

In recent years a significant attention has been paid to experimental and theoretical study of the charge carrier 
transport in disordered organic materials (for good reviews see Refs. [0, El p[). The obtained results profoundly 
changed our understanding of the characteristics of the charge carrier transport in these materials. The essence of a 
modern paradigm of the charge carrier transport is a close connection between spatial correlation of energy levels U{r) 
of transport sites in particular material and its transport properties [0, o] . The mostly well understood reason for a 
strong spatial correlation is a presence in the material some particular kind of molecules, e.g. dipolar or quadrupolar 
molecules y, |^, ^. The most important parameter of the theory is the variance a^ ~ (C/^) — {U) which gives a scale 
of energy fluctuations in the material. For example, in dipolar materials mobility has the following dependence on 
electric field E and temperature T 



H = ^lo exp 



Aw)--' 



/ a n3/2 



eaE 



(1) 



where computer simulation for 3D model gives A ~ 9/25, Cq — 0.78, F = 2 B, while an exact solution for ID case gives 
A = 1, Co = 2, F = Q. The particular kind of mobility field dependence In/i ex E^^'^ (a Poole-Frenkel dependence) 
is a result of slow decay of correlation function C{r) = (C/(r)C/(0)) oc 1/r. This general theory describes well all major 
features of photoinduced charge transport in disordered organic materials. 

Until now a very limited attention has been paid to the influence of charges on the correlation properties of disordered 
organic materials (the only exception is the paper M). These charges may be movable charges flowing through the 
sample or static charges inserted into the sample, e.g. by the pulse of accelerated electrons. Effect of charged particles 
on transport properties was mainly considered using the framework of the oversimplified Poole-Frenkel model of the 
isolated Coulomb trap. This approach totally neglects collective effects which are significant because of the long 
range nature of the Coulomb potential and, thus, is limited to the case of extremely small concentration of charged 
particles. This very long range nature of Coulomb interaction requires to include into consideration the finite thickness 
of transport layer. In the most typical configuration transport layer consists of a film of organic material bounded 
by highly conductive (metallic) electrodes. We show that contribution of image charges created by these electrodes 
significantly modifies properties of organic medium. For example, this effect removes a notorious Coulomb divergence 
and makes possible to calculate in a finite form properties (mean carrier energy, its variance, correlation function, etc) 
of organic medium with nonzero total charge density. Hence, this approach presents a first step in order to understand 
the influence of finite concentration of carriers on the transport properties of materials (beyond the effect of space 
charge) . 

In this paper we limit our consideration exclusively to the case of a single charge carrier moving in the presence of 
static charge background, thus avoiding a tremendous problem of solving a self-consistent problem of the simultaneous 
motion of many charges. Additionally, the problem of static charges presents a significant interest by itself, because 
this very situation may realize in some particular cases, for example in the case of radiation- induced charge generation. 
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It will be shown also that if surface of the conducting electrodes is rough, then surface charge accumulated under 
voltage applied generates an additional disorder in the bulk of organic transport layer. We will consider its influence 
on transport properties of the device. 

At present we have two major approaches to consider charge transport in disordered organic materials. One is a 
computer simulation of the hopping process, and another is an exact solution for the equilibrium carrier mobility in 
a simplified ID model [Q 

IJ'=^S . (2) 

e/d2/exp(-ey + /?2[C7(0)-C(ay)]) 



where e = ea(3E, j3 = 1/kT, and a is a minimal charge-dipole separation. In the most thoroughly studied case - 
charge transport in dipolar glasses - it was shown that these two approaches give essentially the same results with 
not very important differences for some numeric coefficients [pi . 

For this reason we use in this paper the ID model for evaluation of the transport properties of the devices. Word 
"devices" instead of "materials" is used on purpose, because for the extremely long ranged Coulomb forces the 
relevant properties of the transport layer (correlation function etc) depend on the size of the device (on the thickness 
of transport layer) as well as on microscopic properties of the transport medium. In this case a magnitude of the 
disorder in some particular domain of transport layer is not a constant but depends on the location of this domain. 
In other words, charge induced disorder present a case of inhomogeneous energy disorder. In this paper we develop a 
generalization of the ID model for this new case. 

II. CHARGE TRANSPORT IN ID INHOMOGENEOUS SYSTEMS 

Until now very little is known about 3D charge carrier transport in inhomogeneous random systems. At the same 
time, the ID model suggested in Ref. || proved to be a useful tool in revealing the dependence of a carrier mobility on 
electric field and temperature g. For this reason it seems to be a reasonable starting point to study charge carrier 
transport in the framework of the ID model. 

We are going to find a stationary solution Cs{z) for the diffusion of a charge carrier in the presence of external 
electric field E, random potential U{z), a stationary source of carriers located at z = 0, and a perfect sink at z — L 

dc dJ r ^r \ /ON 

^ + ^=V(.), (3) 



J{z,t) ^-D 



dc ^ fdU 



J{0,t) = 0, c{L,t)^0, P=l/kT. 



The most important quantity to calculate is a stationary velocity 

laL 



Jo dzcs[z) 



A direct solution of Eq. (0) gives 

LD 



/o d2e-/3[C/(^)-e£;^] J^ ^^,^(3[U(z')-eEz'] 



(5) 



Unfortunately, to perform the subsequent averaging of v over random field U{z) is a difficult task. For this reason we 
will perform the calculation of the average time {ti) for the carrier to reach the opposite electrode, where t^ = L/v. 
One may reasonably assume that functional dependences of (u) and L/ (i^) on the relevant parameters {E,T, etc) 
should be the same. We have 



rL rL 



{tL) = ^f dx f dzV''^(^-^')(e'^f^(^')-^Wl). (6) 

If U{z) is a Gaussian random field, then 

{tL) = ^l'^dz j^^ dz'e'='^^(^-^') exp (^/3 [{U{z')) - {U{z))] + \(i^ {[U{z') U{z) - {U{z')) + ([/(z))]^)) (7) 



D 



L pL 

dz 



d,'e-m^-^') exp (^(3 [{U{z')) - {U{z))] + \(i^ [C(z, z) + C{z' , z') ~ 2C{z, z')]\ , 



C{z,z') = {U(z)U{z'))~{U{z)){U{z')). 

In the case of L ^ oo and a homogeneous random energy Uiz), when {U{z)) — const and C{z,z') depends on 
z — z' only, this result exactly corresponds to the result of Ref. ^. Hence, result (0) is a direct generalization of the 
corresponding result of Ref. to the case of inhomogeneous random potential. 

Transformation to new coordinates s = z + z',q = z — z' reveals an important symmetry of Eq. (Q) in the case 
when {U{z)) = 



<^-> - Jd 



ds / dqe'"^'^'ig{s, q)+ ds 

lo Jo Jl Jo 



dqe'^'^'^g{s,q) 



(8) 



g{s,q) =exp<{ -(3'^ 



C 



s + q s + q 



C 



s — q s — q 



2C 



s + q s — q 



2 2 J \ 2 ' 2 J V22 

Making an additional transformation s => 2L — s in the second integral in Eq. (^ we have 



<"' ^ ^ (i' 



ds I dqe"'^^'' 





g(s,g)+ / ds f dqe'''P^^g{2L^s,q) 
Jo Jo J 



(9) 



2D 



ds / dqe'^^'^'i [g{s, q) + g{2L - s, -q)h{2L - s, q)] . 



(note that g{s, q) is an even function of q). The sum g{s, q) + g{2L — s, —q) in Eq. (m (and, hence, (ii)) is invariant 
under the transformation z^L~z,z'^L~z', so 



(tL) = {tt) , 



(10) 



where (i^) corresponds to carriers generated at z = and absorbed at 2; = i, and (i^) corresponds to the reverse 
process (with the correspondent inversion of the electric field). We will see that in spite of this symmetry the temporal 
dependence of the photocurrent transient may be very different for these two cases. 

III. GENERAL PROPERTIES OF THE CHARGE-INDUCED ENERGY DISORDER 
A. Potential energy, generated by random sources 



We consider the most simple model of the transport layer - a slab of organic material (infinite in two dimensions) 
bounded at 2; = 0, i by two parallel conducting electrodes having zero potential. Organic medium is modeled by 
the regular lattice with lattice constant a and fraction c of all sites occupied by point charges. Since we are mostly 
interested in the case c ^ 1, so the model of the perfect space lattice is not a significant limitation. Lattice sites 
are also occupied by organic molecules (of dipolar, quadrupolar nature etc), their random positions and orientations 
produce an intrinsic energy disorder, which we neglect in the most part of the paper, because it is independent of the 
charge induced disorder, provides an additive contribution to quantities in question (correlation functions etc), and 
its properties have been extensively studied ||, |, |, 0, |, ||, |l^, [l|, [l2[ ^ ^ . 

Let us consider a general case of a potential energy, generated by random sources 



U{f) = ^ rinf{r, f„, j/„). 



(11) 



Here source function f{r,fn,iyn) describes the contribution from a particular source to the total sum (11) and i^ is 
an additional random parameter (or parameters) that may affect the strength of the contribution (e.g., it may be an 
orientation of dipole in the case of dipolar sources). Note that summation in Eq. (|ll|) is performed over all lattice 



points and rjn — 1 ii a, real source is located at the site n and rjn — otherwise. In this paper we need to calculate 
positional averages of two kinds 

{(VnAil^n))) - iVn) {A{l^n)) = C {A{u^)) , (12) 

and 

{{VmVnA{l'rn)A{iyn))) = (vl) Srnn (^^(t^n)) + (??m) {Vn) (1 " S^n) {A{ly„^)) {A{l^n)) (13) 

= cSran {A"^ (l^n)) + C^(l ^ Smn) (^(i^m)) {A{Vn)) ■ 



Double brackets in the left part of Eqs. (|lj) and (g^ denote a simultaneous averaging over positions and parameter 
v. 

B. Source function and correlation function for point charges 

In the case of random potential energy generated by randomly located point charges function f{r,fn,i^n) is pro- 
portional to the electrostatic potential of the charge located at r = r„ (here and later we assume that charge carrier 
has a positive charge) 

/(r, r„, I'n) = eLp{f, r„), (14) 

and (p{r, fvi) obeys a Poisson equation 

with zero boundary conditions 

V'\.=o^v\.=L=0- (16) 

Hence, ip{r, r„) is proportional to the Green function of the Laplace equation with the corresponding boundary 
conditions 

^(r,r„) = -^G(f,r;). (17) 

Green function G{f, r7i) is calculated in the Appendix H and has a form 

^/^ ^ \ 1 V^ . STTZ . STTZn ^^ /STTflX , „, 

^ ' "^ " "^ ^ T ' ~L^^" \l) ' ^^^ 

where iiro(x) is the Macdonald function and p"^ = {x — XnY + [v ~ VnY ■ 

We consider here a slightly more general case when two types of static charges are present: positive with fraction 
c+ and negative with fraction c^ (we assume that these fractions may vary with n). Note, that c^ +c^ < 1. Average 
energy of the charge carrier is 

{Uir))^-^J2(^+^c-)G{r,f„), (19) 

n 

and the correlation function C{r, r') is 

Cir,P) ^ (f/(r)C/(rO) - (Uir)) (ui?)) = ^^ ^ [c+ + c" - (c+ - 0")^] G(r, f„)G(P, f„). (20) 



C. Uniform distribution of charges 

As a simple example we provide a result for the case of uniform charge distribution with c^ ,c^ — const. Here 

(f/(,)) = ^^£!(£!_£l,(L-z), (21) 



and 



s L L ' \ L J ' 

s=l 



(22) 



s=l 

In the particular case p = 

4716^ [c+ + C" - (C+ - C-)^] 

C{z,z',0) — ^-^5 z_(L — z+), Z-|_ = maxfz, 2'), z_ = min(z, z'). (23) 

Equation (03) shows that even in the case of uniform charge distribution the resulting energetic disorder is not a 
homogeneous one: the correlation function does depend on z, z' separately, and not on the difference z — z' only, as it 
should be for the homogeneous case. Of course, the reason for this phenomenon is an effect of conducting electrodes, 
which is important for a very long range Coulomb interaction. Magnitude of the disorder induced in a transport layer 
could be measured by the variance <j'^{z) = C{z, z). For the uniform distribution of charges 

'^'(^) = ^2^3;^ -^(L - z), (24) 

it attains a maximum at z — L/2 and goes to zero at the electrode surfaces (as it should be, because electrodes have 
a constant potential). 

A direct calculation gives for the time (0) 

(^i) - T7 / d(J di' exp [A?5(C, O] , (25) 



D JO 



(26) 



A = , I = —, f = 2tt\c2 — 2ttIci , ci = c — c^ , C2 = c + c^ — (c — c^ ) . 

ea a ' e 

Introducing new coordinates C — S, + ^' Sind rj — ^' — £^ we have 

(iL> = ^ / dije^], [\l (-2^Ac2r;2 - /y^)] / dCcxp {-2ttI^ Xc^Qt^) 

" 2D^P\ci I ~ """"P [2^A2?c2(r/ - r;^) - ePELr^] sinh [2^PAci {r, - r;^)] . 

One can see that Eq. ( p6| ) is an even function of ci , so charge carrier interacts in the same way with static backgrounds 
created either by trapped electrons or holes. It is easy to show that this particular property is valid for any charge 
distribution having symmetry {U{z)) = {U{L — z)) (spatial symmetry of C(z, z') is not important because it is an 
even function of ci). We should emphasize that this particular property holds for the Gaussian approximation (M) 
only. 

Let us discuss Eq. (Eq) in the particular case Ci = 0, i.e. when the medium is neutral. This case was briefly 
discussed in Ref. for the infinite medium. We have 

(iL)„ - -^ y dr,(l ~ ,7) exp [XI {~2nXc2v' - fv)] - H{f) + j^^' (27) 



Hif) - ^ y dr/exp [a; {-2ttXc2V^ - fv)] 



exp 



flP 



Finally, 



(^i)n 



L^ 



A^DX^lc2 



2DX^/2k^ \8t:c2 



{A + B)e^ [erf (A + B) - erf (B)] 



erf 



/'4:TT\lC2+lf 

V 2y/2Tdc^ 



erf 



f If 



\2^/2ttIc^ 



1 -e 



-/l^-2AB 



(28) 



A = Xx/2ttIc2, B = f 



8tTC2 



The most prominent featm'e of Eq. (pi) is its behavior for L -^ cx3. If / > (i.e. if £■ > Ec,-), then (iL)„ oc i and 
carriers have a well defined average velocity which does not depend on L 



e(3D {E ~ E^r) 



^rr 



2TT(3e^C2 



(29) 



If / < 0, then log(ii)j^ c>c L and average carrier velocity goes to 0. Transition occurs at i? = E^r in full agreement 
with Ref. 0. 

We would like to note a close connection between the case of charge carrier transport in the infinite neutral medium 
c^ and the anomalous diffusion of the particle in a ID random force field, described by the Langevin 



,+ _ 



with c 
equation 



with a thermal noise ri(t) 



dz 
'dt 



-i^(z)+ry(t). 



kT 

(ry(i))^0, {rj{t)rj{t'))^2—S{t~t'), 

7 

and a quenched random force F{z) being a Gaussian white noise 

(F(z)) = ^0, {Fiz)Fiz')) - F^ = a5{z - z'). 



(30) 



(31) 



(32) 



This problem was discussed thoroughly in Refs. nsl E6l it demonstrates the same transition from the mobile carriers 
with non-zero average velocity to the immobile ones. The relation between the case c"*" = c~ (for L -^ oo) and the 
problem ( pO[ ) arises from the identical forms of the corresponding energy-energy correlation functions. In both cases 

( \U(^z) — 'U(z')\ ^ ex |z — z'\. A more careful analysis of this relationship will be studied elsewhere. 

In the general case of macroscopically charged medium with ci 7^ we can estimate time (R) for strong disorder, 
when 27r^^A|ci| ^ 1. In this case we can retain only the leading exponential term in (E9), so 



(iz 



2A/(Ac2-H/|ci|) 



2DI'^X\cy\ V 27rA/(Ac2 + ;|ci|)-e/3£:L 



exp 



{2-nXl{Xc2^^l\cx\)-efiELy 
87rA/(Ac2-f Z|ci|) 



(33) 



This expression is valid if 



27rA/(Ac2 + /|ci|) - e(iEL > [27rA?(Ac2 + ?|ci|)] 



1/2 



Note, that if ci 7^ 0, then log (t^) ex L^ for L — > 00, so for a infinite charged medium there is no non-dispersive 
transport. Field dependence of time ( P6|) is shown in Fig. |l|. 

The case of the uniform distribution of static charges does realize in the case of radiation-induced charge generation 
in disorder organic materials. In this particular case movable charges are induced in the material by the pulse of 
accelerated electrons. In the most typical case electrons have high energy, so concentration of electron-hole pairs 
(usually, electrons get trapped and holes are movable) created in the bulk of the material is approximately constant 
and does not decay while going away from electrode. We could expect that the model of uniformly charged material 
is a reasonable model for description of transport properties in the radiation-induced generation, at least for large 
times, when concentration of movable charges is much smaller than the concentration of static charges. 

A different case of static charges trapped at the close vicinity of the electrode (this case does realize for the 
radiation-induced generation of charge carriers by the pulse of low energy electrons) will be studied elsewhere. 
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FIG. 1: Field dependence of (t^) for the case of uniformly charged background with L/a = 150, e^jijea = 20, 
different values ofc^: lxl0~^, 3x10^"", and 5 x 10^^, from upper curve downwards, correspondingly. 
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D. How valid is the Gaussian approximation? 

Particular case of the uniform charge distribution gives us a good possibihty to discuss a vahdity of the Gaussian 
approximation for the ID transport model. Gaussian approximation is valid for the main body of the distribution 
if (7^(z) ^ Ul^ where C/g = e^c^l'^ jea is a typical contribution from the single charge source, thus ensuring that the 
resulting potential energy is not dominating by the single contribution. This inequality leads to 



47rci/3 
ah 



z(L- z) > 1. 



(34) 



Hence, the Gaussian approximation is valid in the bulk of the layer if irc^i^L/a ^ 1, but in the near vicinity of 
electrode at z < Zb, Zb ~ a/^.iTc^/^ non-Gaussian effects are significant. If 'iirc^/^ > 1 (c > 5 x 10""*), then for the 
whole transport layer the main body of distribution has a Gaussian shape. 

Nevertheless, a Gaussian shape for the tail of the distribution holds only if a more strong inequality is valid: 



\z) » U^ 



where C/„ 



/ea is a maximal contribution from the single source. This leads to 



Attc 
aL 



z(L- z) > 1, 



(35) 



so tail of the distribution is Gaussian in the bulk of the layer if -KcL/a 3> 1, and at the vicinity of electrodes if Attc > 1 
(or c > 0.1). The last condition seems to be difhcult to match. 

IV. RANDOM POTENTIAL ENERGY, GENERATED BY CHARGED ROUGH METALLIC SURFACE 

A. Basic formalism 



A typical transport device can be described as a flat capacitor with the region between conducting electrodes filled 
by an organic material (transport layer with the thickness L). If voltage Vq is applied to the capacitor, then surfaces 
of the electrodes develop a surface charge with the density Q proportional to the electric field E = —Vq/L, namely 
Q — E/Att. In the ideal case of absolutely flat electrodes' surfaces the surface charge just maintains this very field 



8 

E = const in the bulk of organic material, according to the well known formula for a correspondent distribution of 
the electrostatic potential ip in the bulk of the capacitor 



(p{f) = 



QdS 
R ' 



(36) 



where integration is carried out over the electrode surfaces and i? is a distance between point r and a particular part 
of the surface. Situation is different in the case of a rough surface of the electrodes (Fig. ||a). In this case we can 
describe surface of the electrode as a random surface, so the distance R in Eq. ( p6| ) fluctuates, thus the resulting 
potential (p is a. random function. Hence, a roughness of the electrode surface induces an additional energetic disorder 
in the bulk of the organic transport layer. We are going to estimate a magnitude of this kind of disorder and its effect 
on the transport properties of devices. One could immediately note two important properties of this kind of disorder: 
its magnitude increases with the increase of the electric field E applied to the device, and decreases while going away 
from the electrodes (so this is another case of inhomogeneous disorder) . 

Let us calculate the potential distribution inside a capacitor having rough metallic electrodes. Surface of one 
electrode obeys an equation z = ho{x,y), and another one z — L + hL{x,y), where ho{x,y) and hL{x,y) are random 
functions and /iq, h^ <C L. We define the locations of the average electrode planes in such a way that {Iii^^l) — 0. We 
are going to study variation of the potential for z ^ h. 

Inside the capacitor we should solve a Laplace equation 



d'^ip d'^Lp d'^ip 
dx^ dy^ '^^^ 



dz^ 



with potential (p obeying the boundary conditions 



V>\z=ho{x,y) - 0' flz^L+hLix.y) " ^0- 



(37) 



(38) 



We are going to find a correction for the potential distribution inside a capacitor which is resulted from the roughness 
of the electrode surfaces 

ip{f)^-Ez + 6ip{f). (39) 

Calculation of the first order correction to the potential distribution with Sip (x h is performed in Appendix ^. 
The energy-energy correlation function has a form (note that {6^{f)) — 0) 

Cir{,ri) = e" {5ip{r{)5ip{r'2)) (40) 



47r2 J sinh^ kL 



UQfto{k) sinhfc(i — zi) sinhfc(i — Z2) + UiilL{k) sinhfczi sinhfcz2 



Uj=eE{h]) 



2\l/2 



%(Pi -P2) = 



>^) 



{h,{p,)h,{p2)) , %(fc) = / dpe-^'^P^M^ J = 0, L 



(here k and p — {x,y) are 2D vectrors). We can simplify equation (Eof) in a situation when the surface correlation 

function i^{p} has a single length scale / (surface correlation length). In a very common situation / ^ {^/i^) (this 
is so, for example, for the so called naturally rough metal surfaces |l^, |l^). li z <^ L, then significant domain in 
integral (^) is fc^ ^ 1. li zi,Z2,l <t^ L (we consider a vicinity of one electrode), then we can simplify equation (^) 



C^(rl,r5) 



f/2 

47r2 



dkflo{k) exp{ik{p{ - p^) - k{zi + Z2)), 



so if zi, Z2 3> ^ then 



C(rl,r^) ccUqI / dkexp{ik{p{ - p'2) - k{zi + Z2)) 



_CWfl+f2) 

[(zi+Z2)2+p- 



13/2 ■ 



(41) 



(42) 



This universal behavior takes place when a large area S ^ l'^ oi the electrode gives a contribution to the distribution 
of the potential. Such regime is observed for zi, Z2 » I only, because essential area S (x z"^. Note, that in this very 
case we should expect that the distribution of Sip has approximately a Gaussian form (because many independent 
areas of the surface give contributions to the resulting distribution). In the opposite case, when zi,Z2 ^ I, then 
C(rl,rl)«C/2. 



B. What should we expect to observe in experiments? 

We can calculate time (t^) for the case of rough electrodes using Eq. (0). For the most interesting case / » (h?) 
results are shown in Fig. 0. The most important feature of the mobility dependences is an abrupt decrease of the 
carrier velocity when E reaches a critical value. Using a saddle-point method one can show that -Ecrit oc l/e(3(h?j. 

To the best of our knowledge, such abrupt drop of the mobility has not been observed in experiments. Quite 
possibly, i^crit is so high that it cannot be reached because of the dielectric breakdown of the layer. If ^ = 15 nm, 

(/i^) — 1.8 nm (data from Ref. p7|), then i^crit ~ 2 x 10^ V/cm. This field is indeed high, but it could be reached 
for layers of high quality (see Ref. |19|). Another possible reason of the discrepancy could be a failure of the ID 
model in this particular case. Finally, our approach totally neglects a correlation between rate of the charge carrier 
injection from some electrode's area and its profile. We may expect that our approach is valid for rather "smooth" 
rough surfaces, where we cannot expect a significant variation of the injection rate over electrode surface. Anyway, 
we believe that a careful study of charge transport in high quality layers with simultaneous control (by the use of 
scanning tunneling microscope) of the roughness of the electrodes could confirm (or disprove) our results. 

We should expect that a significant effect of the roughness could be observed for weaker electric field if we consider 
a temporal behavior of the photocurrent lit) for the time-of-flight experiments. At present we cannot produce any 
analytic or simulation results for I{t). Anyway, we may safely suppose that: 

1) Transport should be more dispersive in the case of rough electrodes in comparison to the case of smooth electrodes, 
just because of the increase of disorder. Additionally, this increase should lead to the decrease of mobility. 

2) If a rough electrode serves as an injector, we should expect an unusual transformation of the temporal behavior 
of the photocurrent I{t) with the increase of the electric field E. Typical distribution of sites' energies in this case 
is shown in Fig. 0b. In a weak field region we should expect a transient of the usual kind (Fig. He). In strong field 
region, with the increase of the magnitude of the disorder near the injecting electrode, transient should transform to 
the form shown in Fig. ||d. A reason for the development of the region where current increases with t is carriers' 
acceleration when they move to the region with smaller disorder. This kind of photocurrent transformation with the 
increase of E was indeed observed in Ref. ^ 

3) We noted previously a principal symmetry of the equation (M), so (^7*) — (^IT) even when electrodes have very 
different degree of roughness (compare the left and right electrode in Fig. ^) . But if we consider a temporal behavior 
of the current I{t), it could be very different for these two cases. For not so small E, if a smooth electrode serves as 
injector, we should expect a usual form of /(f) (see Fig. gc), but if a rough electrode serves as injector, we should 
expect the form shown in Fig. pp. 

V. CONCLUSIONS 

We suggested a new approach for calculation of statistical properties of disordered organic materials with nonzero 
charge density. Explicit consideration of the finite transport layer bounded by conducting electrodes removes a long 
range Coulomb divergence and gives finite values for the relevant characteristics of the device (correlation function 
etc). We presented a generalization of the ID model of the charge carrier transport in disordered organic matrices 
to the case of inhomogeneous disorder. This approach permitted us to calculate transport properties of matrices 
with charge induced disorder. A similar approach permits us to calculate an energy-energy correlation function for 
the transport layer bounded by rough conducting electrodes. We expect that the roughness of the electrode surface 
should lead to the unusual temporal behavior of the photocurrent transient (appearance of the region where current 
increases with time) in the case of strong electric field. 
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APPENDIX A: CALCULATION OF THE GREEN FUNCTION FOR THE LAPLACE EQUATION 

Taking into account a geometry of the problem, we may write 

G{r, ^") = 4^ / dke'^^f^Gt^i^z, z„), (Al) 
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where k and p — {x — x„, y — j/„) are 2D vectors, and function Gk{z, z„) obeys an equation 

'^^'^ PGk = 6{z-Zn), Gfe(0,z„)=Gfc(L,z„) = 0, Gfe(z,z„) = Gfe(z„,z). (A2) 



dz^ 
A direct calculation gives 

sinhfc2;_sinhfc(L - z+) • / n / a o^ 

Gfc(z,z„) = , . ■ , ^ , z+ = max(z,z„), z_ = min(z, z„). (A3) 

We would like to note another useful expression for Gk{z, Zn) which is based on the series of the eigenf unctions of 
the operator d'^ — k^ 






From Eq. (A4) it immediately follows that 



'P{r, r„) = ^ 2^ sm — sm —j—K^ VJT) ' ^ ' 



s=l 

where Kq(x) is the Macdonald function. Last expression is useful to show behavior of (y9(r, r„) for p^ L. Presence 
of metallic electrodes transforms the usual Coulomb decay to the exponential one 

^ ^ 4e nz nzn L f np\ 

<^('',?'n) ~ ^sm — sm— W — exp(^-— J , p > L, (A6) 

thus removing a nasty long range Coulomb divergence. 

APPENDIX B: FIRST ORDER CORRECTION TO THE POTENTIAL DISTRIBUTION FOR A 
TRANSPORT LAYER BOUNDED BY ROUGH ELECTRODES 

To solve Eq. ( |37| ) let us make a transformation of variables 

X = x, Y = y, Z^J^j^-^-, (Bl) 

SO the boundary conditions ( |38[) transform to 

^\z=o = 0' ^lz=L = Vo. (B2) 

In new variables {X, Y, Z) the Laplace equation (^^ has a form 

f d^Lp dZ , 9V dZ\ , dip fd^Z , d'^Z" 

(B3) 



9^</j 9^(/3 9^(/3 


'fdzv /dzy /dzy 

\dxj \dy) \dz) _ 



"'^'"'^'"' ^^ (yaxaz 9a; ^ ar9z 9y j ^ az (^ 9x2 + ay2 j o 



We will seek a solution of Eq. (B3) by the perturbation theory (a formal small parameter being h) 

(^ = ^(/3„, (^„~0(/i"), ipa^VaZ/L, •^n\z=Q^ Vn\z=L ^ ^-^ n > 1. (B4) 

n 

Our main goal is to calculate a leading term for the correlation function C{f{, f2) — ^ {S'p{fi)dip{T^)) . For this reason 
we should calculate ipi only. 

Equation for an every term tpn has a general structure 

Aipn = Jn, (B5) 
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where the source J„ depends oniy9fc,0<fc<n — 1. Solution of Eq. (B5) is 

Vn {r) = / df{G{r, f{ ) J„ (rl ) , 



here G{rTfi) is the Green function ( |A1[ ) for the Laplace operator with zero boundary conditions (B4) 
The source term for the first order correction is 



Ji = 



L 



z\ z 



A, = 



d^ 



d^ 



dX^ dY^' 



(B6) 



(B7) 



ipi{f) = -r I di\G{r,f{) 



Calculation of the integral (B8) gives 



"Piir) 



Att^L 



1 - ^^ A^hoipl) + ^A^hLipi) 



dk / dpi / dZie'''^P-P''h'^Gk{Z,Zi) 



1 - §) hoipi) + ^hL{p{) 



(B8) 



(B9) 



Vn 



Att^L 



dke^'^P dZik^Gk{Z,Zi) 



l-j^)hS) + ^hL[k) 



47r2L 



dke'^P 



Z sinh fc(L - Z) \ -. fZ sinhfcZ 
L sinh kL J \ L sinh kL 



hL{k) 



The total first order correction to the potential distribution (here we are going back to the original coordinates) is 

<5^ - (^ - ((^) = --^ [(l - I) h^{x, y) + ^hL{x, y)] (BIO) 






dke'^P 



l_ L - sinhfc(L-z) '. ^ ^ 
L sinh fcL 



z sinh kz 



L sinh fcL 



/iL(fc) 



'Att^L 



dke^'^P 



sinhfc(i — z] 
sinh kL 



hoik) 



sinhfcz 
sinh kL 



hdk) 



[1 

[2; 
[3; 

[6; 

[9 

[lo; 
[11 

[12 
[13 
[14 
[15 
[16 

[17 
[18 
[19 
[20 



P. M. Borsenberger and D. S. Weiss, Organic Photoreceptors for Imaging Systems (Marcel Dekker, New York, 1997). 

H. Bassler, Phys. Status Solidi B 175, 15 (1993). 

M. Pope and C. E. Swenberg, Electronic Processes m Organic Crystals and Polymers (Oxford University Press, New York, 

1999). 

D. H. Dunlap, P. E. Parris, and V. M. Kenkre, Phys. Rev. Lett. 77, 542 (1996). 

S. V. Novikov, D. H. Dunlap, V. M. Kenkre, P. E. Parris, and A. V. Vannikov, Phys. Rev. Lett. 81, 4472 (1998). 

S. V. Novikov, D. H. Dunlap, and V. M. Kenkre, Proc. SPIE 3471, 181 (1998). 

D. H. Dunlap and S. V. Novikov, Proc. SPIE 3144, 80 (1997). 

A. Dieckmann, H. Bassler, and P.M. Borsenberger, J. Chem. Phys. 99, 8136 (1993). 

S. V. Novikov and A. V. Vannikov, Russ. JETP 79, 482 (1994). 

R. H. Young, Philos. Mag. B 72, 435 (1995). 

S. V. Novikov and A. V. Vannikov, J. Phys. Chem. 99, 14573 (1995). 

S. V. Novikov and A. V. Vannikov, Proc. SPIE 2850, 130 (1996). 

S. V. Novikov and A. V. Vannikov, Synth. Met. 85, 1167 (1997). 

S. V. Novikov and A.V. Vannikov, Proc. SPIE 3144, 100 (1997). 

J. P. Bouchaud and A. Geogres, Phys. Rep. 195, 127 (1990). 

J. P. Bouchaud, A. Comtet, A. Geogres, and P. Le Doussal, Ann. Phys. 201, 285 (1990). 

M. Rasigni, G. Rasigni, J. P. Palmari, and A. Lebaria, Phys. Rev. B 23, 527 (1981). 

G. Rasigni, F. Varnier, M. Rasigni, et al, Phys. Rev. B 27, 819 (1983). 

L. B. Schein, A. Peled, and D. Glatz, J. Appl. Phys. 66, 686 (1989). 

M. Sugiuchi and H. Nishizawa, J. Imaging Sci. 37, 245 (1993). 



12 







■^tl^ ^ 










♦ ♦ 




* , ♦ 


♦ ♦ 






♦ ♦ ♦ ♦ ♦ ♦ ♦♦♦ 


3 


♦ ♦ « 

♦ 

* 


♦ 


* 


< <t2> 




• 






♦ 





(a) 



(b) 



*♦♦♦♦♦..♦♦♦♦•*♦....*♦♦♦**♦. 



'♦♦♦.♦. 



■ 




♦ 


o 






♦ 




♦♦♦♦♦♦*♦ 

.♦♦♦♦♦ 

. .♦♦♦♦♦ 
♦♦♦*♦♦ ♦ 




♦ 




♦ 




*♦♦♦♦.. 



time 
(c) 



time 
(d) 



FIG. 2: (a) Transport layer bounded by rough electrodes (note the different roughness of the left and right electrode), (b) A 
typical distribution of the random energy U{z) in the bulk of a transport layer induced by the rough surface of the electrodes in 
the case of significant difference in their roughness. According to Eq. (Iicl) (ti) — (^2), but time dependence of the photocurrent 
may be quite different for these two cases. (c,d) Typical time dependence of the photocurrent in the time-of-flight experiment 
when a smooth electrode serves as an injector (c), or a rough one serves as an injector (d). 
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FIG. 3: Field dependence of (iz,) for L/a = 100, {hlf^'^ /a = 5, {hiy^ = 



bounded by rough metallic electrodes with different values of lo/a: 
Surface correlation function was chosen in the form Qo{p) = exp 



[0. If a = 1 nm and a^ = 0.1 eV, then eaE/aq Ri 1 for £ = 1 x lO'' V/c 



0, in the case of quadrupolar glass with aqf3 = 2 
50, 25, 10, 5, from the right curve to the left, correspondingly. 



^/2lo), which is suitable for a naturally rough metal surface 



